Competition between paramagnetism and diamagnetism in charged Fermi gases 
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The charged Fermi gas with a small Lande-factor g is expected to be diamagnetic, while that with 
a larger g could be paramagnetic. We calculate the critical value of the (j-factor which separates 
the dia- and para-magnetic regions. In the weak-field limit, g^ has the same value both at high and 
low temperatures, gc = l/y/vl. Nevertheless, gc increases with the temperature reducing in finite 
magnetic fields. We also compare the gc value of Fermi gases with those of Boltzmann and Bose 
gases, supposing the particle has three Zeeman levels a — ±1, 0, and find that gc of Bose and Fermi 
gases is larger and smaller than that of Boltzmann gases, respectively. 
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I. INTRODUCTION 

Magnetism of electron gases has been considerably 
studied in condensed matter physics. In magnetic field, 
the magnetization of a free electron gas consists of two 
independent parts. The spin magnetic moment of elec- 
trons results in the paramagnetic part (the Pauli para- 
magnetism), while the orbital motion due to charge de- 
gree of freedom in magnetic field induces the diamagnetic 
part (the Landau diamagnetism) P] . The Pauli paramag- 
netism and the Landau diamagnetism compete with each 
other. For electrons whose Lande-factor g = 2, the zero- 
field paramagnetic susceptibility is two times stronger 
than the diamagnetic susceptibility, so the free electron 
gas exhibits paramagnetism altogether. 

Magnetic properties of relativistic Fermi gases have 
also been under extensive investigation. Daicic et al. 
developed statistical mechanics for the magnetized pair- 
fermion gases and found that the intrinsic spin causes 
important effects upon the relativistic para- and dia- 
magnetism 

The study of ultracold atoms has stimulated renewed 
research interest in the magnetism of quantum gaseslsl- 
0). When atomic gases are confined in optical traps 
their spin degree of freedom becomes active, leading 
to the manifestation of magnetism. Theoretically, the 
paramagnetism [3] and ferromagnetism (8l-[l0| in a neutral 
spin-1 Bose gas have ever been studied. In experiments, 
magnetic domains have been directly observed in *^Rb 
condensate, a typical ferromagnetic spinor condensate [5]. 
Very recently, the exploration of magnetism in quantum 
gases has been extended to the Fermi gas. It is already 
observed that an ultracold two-component Fermi gas of 
neutral ^Li atoms exhibits ferromagnetism caused by re- 
pulsive interactions between atoms l6|. 

Furthermore, it is possible to realize charged quantum 
gases from neutral ultracold atoms. So far, cold plasma 
has been created by photoionization of cold atoms [TT|. 
The temperatures of electrons and ions are as low as 100 



mK and 10 /iK, respectively. The ions can be regarded 
as charged Bose or Fermi gases. Once the quantum gas 
has both the spin and charge degrees of freedom, there 
arises the competition between paramagnetism and dia- 
magnetism, as in electrons. Different from electrons, the 
g-factor for different magnetic ions is diverse, ranging 
from to 2. 

The y-factor is a characteristic parameter to measure 
the intensity of paramagnetic effect. It is expected that 
the quantum gas displays diamagnetism in small g region, 
but paramagnetism in large g region. The main purpose 
of this paper is to calculate the critical value of g. We 
also present a comparison with results of charge d sp in-1 
Bose gases which have been obtained previously [l2|- 

II. THE MODEL 

We consider a charged Fermi gas of N particles, with 
the effective Hamiltonian written as 

H - fiN = Dl ^ (cjfe, + - /i) n-jk.a, (1) 

where is the chemical potential and £jk^ is the quan- 
tized Landau energy in the magnetic field B, 
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with j = 0, 1, 2, . . . labeling different Landau levels and 
oj = qB/{m*c) being the gyromagnetic frequency of a 
fermion with charge q and mass m* . = qBSxy/ (27r?ic) 
marks the degeneracy of each Landau level with Sxy the 
total section area in x, y directions of the system, 
denotes the Zeeman energy, 
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where g is the Lande-factor and a refers to the spin-z 
index of Zeeman state 1^", cr). 

The grand thermodynamic potential of the Fermi gas 
is expressed as 

^T^o = -^DlY^ ln[l + e-'5(^^-^+^'-^)], (4) 
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where P = (ksT) ^. Performing Taylor expansions and 
then integrating out , Eq. becomes 
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where V is the volume of the system. For simplicity, the 
following notation is introduced, 



F:^[a,S]=Y. 
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where x = fShu and r]„ 
Eq. ^ is rewritten as 



{nuj/2- fi + e„)/{hLu). Then 
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For a system with the given particle density n — N/V, 
the chemical potential is obtained according to the fol- 
lowing equation 
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A similar treatment has been employed to study diamag- 
netism of the charged spinless Bose gas[13| and extended 
to the study of competition of diamagnetism and param- 
agnetism in charged spin-1 Bose gases This method 
is more applicable at relatively high temperatures. 

To determine gc, we need calculate the magnetization 
density M as a function of the magnetic field B and tem- 
perature T. The system is paramagnetic when M > 
while diamagnetic when M < 0. M can be derived from 
the thermodynamic potential by the standard procedure 
Mt^o = (IS)T,y' ^hich yields 



m*c \2T:l3ff 
1. 
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For carrying out numerical calculations, it is useful 
to make the parameters dimensionless, such as M = 
m*cM/{nhq), uJ = fuv/{kBT*), and t = T/T*. Here T* 
is the characteristic temperature of the system, which is 
given by fcsT* = 27rS^n3 /m*. Then we have x = uJ/t, 
rju = 1/2 — ^' /Zd— ga and /i' — ^i/{kBT*) in the notation 
F^[a,S\. Equations ([8|) and ([9]) can be re-expressed as 



1 = cjti/2^i^f [-1,0 



(10) 
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FIG. 1: (a) The total magnetization density (Af), (b) the 
paramagnetization density {Mp), and (c) the diamagnetiza- 
tion density {Md) as a function oi g &t t — 0.2. The dotted 
and dashed lines correspond to Zj = 5 and 10, respectively. 



and 



Mt^o =t3/2 ^(f^ [-3, 0] + x{ga - \] 



F^[-l,Q\-xF^[-l,l] 
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In following discussions, uj is used to indicate the magni- 
tude of magnetic field since it is proportional to B. 



III. RESULTS AND DISCUSSIONS 

First, we look at a spin-^ Fermi gas, setting cr = ±1 
to present the two Zeeman levels. The dimensionless 
magnetization density M as a function of g is shown in 
Fig. [ija). As expected, M is negative in the small g 
region, which means that the diamagnetism dominates. 
For each given value of w, M grows with g and changes 
its sign from negative to positive in the larger g region, 
indicating that the paramagnetic effect is enhanced due 
to increase oig. This phenomenon is also observed in the 
Bose system [12|. gc is just the value of g where Af = 0. 
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FIG. 2: Plots of the critical value of g-factor as a function of 
1/t. The magnetic field is chosen as Zj = 20 (solid line), 10 
(dashed line) and 5 (dotted line), respectively. 



The total magnetization M shown in Fig. [Ija) con- 
sists of both the paramagnetic and diamagnetic contri- 
butions. Figure [ijb) depicts the pure paramagnetic con- 
tribution to M, which is calculated by Mp = gm where 
771 = Til — Figure [ijc) plots the diamagnetic con- 

tribution to M, Md = M — Mp. For each fixed value 
of UJ, the diamagnetization is slightly weakened with 
increasing g. Interestingly, for a charged spin-1 Bose gas 
in a constant magnetic field, the diamagnetism is not 
suppressed but enhanced as g becomes larger [l2j. Com- 
paring Figs[TJa), [Ifb) and [Tfc), it can be seen clearly 
that the increase of M with g is mainly owing to the 
paramagnetization Mp. 

gc is an important parameter to describe the competi- 
tion between the diamagnetism and paramagnetism. Fig- 
ure [5] plots gc for the charged spin-i Fermi gas. The tem- 
perature is described by 1/t and the magnetic fields are 
chosen to be relatively larger, since the method adopted 
here is more applicable at higher temperatures or in 
stronger fields. In the high temperature limit, gc tends 
to a universal value, gclt-s-oo = 0.28868, regardless of the 
strength of the magnetic field. In a fixed magnetic field, 
gc increases monotonically as the temperature falls down. 
But the trend of increasing is slowed down if the magnetic 
field is weakened. 

Present method can not produce valid results in the 
low temperature region if a small magnetic field is chosen. 
In this case that guJ ^ t <C the grand thermodynamic 
potential in Eq. (j4]) can be calculated using the Euler- 
Maclaurin formula 



T.^ij + ^)^ Hx)dx + -i:'{0), (12) 
j=0 •'° 

and then Eq. (|4]) is transformed into 
"t#o = -^E/|W + S(M^E/^W' (13) 
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FIG. 3: The Qc — 1/f curves for charged spin-1 gases obey- 
ing the Bose-Einstein (BE, dotted line). Maxwell- Boltzmann 
(MB, solid line) and Fermi-Dirac (FD, dashed line) statistics, 
respectively. The magnetic field is chosen as Zj = 10 and 5. 



where A = hj3^ / {2'n:m*Y^ , z — e^^^ ^"^ and the Fermi- 
Dirac integral fn{z) is normally defined as 



fn{z) 
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where r(?7,) is a usual gamma function and x = /3e. Then 
the magnetization density, M, is obtained from grand 
thermodynamic potential in Eq. ([Ti 
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where z' = e'^^ +9°''^)/*. After some algebra, we get from 
Eq. dm that 5c|zzj->o = 1/VT2 ~ 0.28868. So in the weak 
field limit, gc has the same value both at the high and 
low temperature limit. This implies that the gc ~ 1/i 
curve is likely to flatten out in the weak-fleld limit. 

To proceed, we make a comparison between Fermi 
and Bose gases. Considering that gc for charged spin- 
1 Bose gas has already been studiedp^. we discuss a 
Fermi gas with three sublevels, a = ±1,0 (a pseudo- 
spin- 1 Fermi gas) to ensure the comparability. Results 
for the spin-1 Boltzmann gas are also obtained. Fig- 
ure [3] shows the gc — 1/t curves for the three kinds of 
charged spin-1 gases. In a given value of w, gc of all the 
three gases displays similar temperature-dependence: gc 
increases monotonously as t reduces. In the high tem- 
perature limit, gc goes to the same value in all magnetic 
fields, gc\t^oo = l/\/8, reflecting that the Bose-Einstein 
(BE) and Fermi-Dirac (FD) statistics coincide with the 
Maxwell-Boltzmann (MB) statistics in this case. 

Figure [3] also demonstrates the difference among the 
three kinds of statistics. For each given fixed value of ZU, 
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the gc—i/t curves of Bose and Fermi gases always locate 
at the two sides of that of the Boltzmann gas. Given the 
same temperature and magnetic field, gc of Fermi gas 
is the smallest. According to our previous research [1^. 
gdt^o = 1/2 for the Boltzmann gas regardless of the 
magnetic field. This means that gc of the Fermi gas does 
never exceed 1/2 in the low temperature, no matter how 
strong the field is. 

IV. SUMMARY 

In summary, we study the interplay between paramag- 
netism and diamagnetism of charged Fermi gases suppos- 
ing that the Lande-factor g is a variable. The gas under- 
goes a shift from diamagnetism to paramagnetism at the 



critical value of g and gc increases monotonically as the 
temperature t decreases in a fixed magnetic field w, and 
the rise in gc is lowered as uJ is reduced. We conjecture 
that gc holds a constant at all temperatures in the weak 
field limit. For a spin-i Fermi gas, gdzj^o = l/\/T2. We 
also briefly compare gc of charged spin-1 gases obeying 
the Fermi-Dirac, Bose-Einstein and Maxwell-Boltzmann 
statistics. The gc — i/t curves of Boltzmann gases are al- 
ways between those of Bose and Fermi gases in the same 
magnetic field. In the high temperature limit, gc of all 
the three gases tends to the same value. 

This work was supported by the Fok Ying Tung Educa- 
tion Foundation of China (No. 101008), the Key Project 
of the Chinese Ministry of Education (No. 109011), and 
the Fundamental Research Funds for the Central Univer- 
sities of China. 



[1] L.D. Landau, E.M. Lifshitz, Statistical Physics. Part 1, 
Butterworth-Heinemann, Oxford, 1980. 

[2] J. Daicic , N.E. Frankel, R.M. Gailis, V. Kowalenko, 
Phys. Rep. 63 (1994) 237 and references therein. 

[3] D.M. Stamper-Kurn, M.R. Andrews, A.P. Chikkatur, S. 
Inouye, H.-J. Miesner, J. Stenger, W. Ketterle, Phys. 
Rev. Lett. 80 (1998) 2027. 

[4] J. Stenger, S. Inouye, D.M. Stamper-Kurn, H.-J. Mies- 
ner, A.P. Chikkatur, W. Ketterle, Nature 396 (1998) 345. 

[5] L.E. Sadler, J.M. Higbie, S.R. Leslie, M. Vengalattore, 
D.M. Stamper-Kurn, Nature 443 (2006) 312. 

[6] G.B. Jo, Y.R. Lee, J.H. Choi, C.A. Christensen, T.H. 
Kim, J.H. Thywissen, D.E. Pritchard, W. Ketterle, Sci- 
ence 325 (2009) 1521. 

[7] K. Yamada, Prog. Theor. Phys. 76 (1982) 443; M.V. 



Simkin, E.G.D. Cohen, Phys. Rev. A 59 (1999) 1528. 
[8] T.-L. Ho, Phys. Rev. Lett. 81 (1998) 742; T. Ohmi, K. 

Machida, J. Phys. Soc. Jpn. 67 (1998) 1822. 
[9] Q. Gu, R.A. Klemm, Phys. Rev. A 68 (2003) 031604(R); 

C. Tao, P. Wang, J. Qin, Q. Gu, Phys. Rev. B 78 (2008) 

134403. 

[10] K. Kis-Szabo, P. Szepfalusy, G. Szirmai, Phys. Rev. A 
72 (2005) 023617; S. Ashhab, J. Low Tempt. Phys. 140 
(2005) 51. 

[11] T.C. Killian, S. Kulin, S.D. Bergeson, L.A. Orozco, C. 

Orzel, S.L. Rolston, Phys. Rev. Lett. 83 (1999) 4776. 
[12] X. Jian, J. Qin, Q. Gu, unpublished. 
[13] G.B. Standen, D.J. Toms, Phys. Rev. E 60 (1999) 5275. 



